Basic Assumptions About Convection 





The calculation of convection requires the we adopt some assump- 
tions about the process. In general, these assumptions are false 
(sometimes by several orders of magnitude), but we are stuck with 
them. 


e There exist discrete “blobs” of material inside a star that have 
sizes of the order of lm, where lm is the mixing length. These 
blobs can move distances of the order of lm before losing their 
identity. (This variable lm is somewhat of a fudge-factor, and many 
of the uncertainties will get hidden in it.) 


e The mixing length is small, compared to other scale lengths in 
the star. In particular, if Ap is the pressure scale height, then 
ln KAP: 





e The blob maintains pressure balance with its surroundings. Thus 
the time for a pressure wave to cross the blob is much shorter than 
the timescale for blob motion, Im/us < lm/v. 


e Acoustic waves, shocks, magnetic fields, rotation, are all negligi- 
ble. 


e The temperatures and densities inside a blob are only slightly 
different than those surrounding it. 


These assumptions are sometimes called “Boussinesq”; they are 
valid for almost incompressible fluids with small density and tem- 
perature variations. 


Convective Energy Transport 


In a radiative region of a star, the temperature gradient is given by 
Vrad; in the limiting case of a fully convective region (i.e., where 
all the energy is transported by convection), the gradient will be 
Vaa. In the outer envelopes of stars, however, the temperature 
gradient can be Vraa > V > Vad. 








The computation of V in the presence of convection requires some 
work. Recall that the flux transported by radiative (and conduc- 
tive) diffusion is 
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When combined with the hydrostatic equilibrium equation (2.2.2) 
and the definition of V (2.4.1), this becomes 
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This gives a relation between Faq and V. Unfortunately, neither 
is known, since Faq carries only part of the star’s flux. The other 
part is carried through blobs of matter, which move with velocity 
v and have density p. If these blobs are in pressure equilibrium 
with their surroundings, then they will be capable of depositing 


Fony ~ pv cp AT (8.1.2) 





where AT is the temperature difference between the blob and the 
surrounding material. Since V+aq is defined as the value needed to 
transport all of the star’s flux, we have 
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The task is therefore to estimate p, v, and AT for the blobs. 
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First, let’s find an expression for AT in terms of the temperature 
gradient V. To do this, consider the behavior of a single blob of 
material. Assume that a typical blob will move a radial distance 
lm before losing its identity, and that, at a given location r inside 
the star, a typical blob will have already moved half this distance. 
At this time, the temperature difference between this typical blob 
and its surroundings will be 
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where the subscripts 2 and s refer to the material inside and outside 
the blob. (We will now drop the subscript s, since those are just 
the global variables of the star.) If we divide each side by T and 
multiply by the pressure scale height, Ap = —dr/dln P, then 
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Note that at this point, the variables V and V; are still unknown. 
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Now let’s estimate the velocity of the blob. According to the defi- 
nition of a, ô, and y, 
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However, if the blob remains in pressure equilibrium with its sur- 
roundings, dP = 0, and if the chemical composition does not 
change, du = 0. Hence 
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which implies 
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This difference in density leads to a net buoyancy force of fẹ = 
—g (pi — p); the work performed by this force while moving the 
blob is then 
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Now suppose that ~ 1/2 this work goes into the kinetic energy of 
the blob while the other half operates on the surroundings. The 
mean velocity of the material passing the blob is then given by 
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With this expression, we can write the convective flux as a function 


of the temperature gradients that are internal and external to the 
blob: 
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Finally, let’s estimate V; for a blob with volume V, surface area 
S, and total heat content (over its surroundings) Q. Clearly, if the 
blob remains in pressure equilibrium, then 


Q= pV cp AT (8.1.8) 


As the blob moves, some of this heat will be radiated away. From 
(3.1.4), the rate at which this will happen will be 
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where x is the direction normal to the surface of the blob. To 
evaluate this, let’s make the simpliest approximation that the tem- 
perature gradient normal to the blob depends on the size of the 
blob, i.e., 
dr T,-T 
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Through (8.1.4), this implies 
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Now, as the blob moves, its internal temperature will change, both 
from the radiative loss of heat, and from adiabatic expansion (or 
contraction). The change in temperature is 
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Now let’s choose some shape for the blob. The ratio of surface 
area to volume for a sphere of diameter lm is S/V = 6/lm, but 
other shapes will have different values. For the moment, let’s just 


assign S/V = Ç/lm. If we now substitute this in, and multiply by 
Ap/T = —(dr/dln P)/T, we get 
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or, if we let ¢ = 9 for compatibility with the literature 
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This last equation provides us the information we need to solve the 
problem. Aside from the local variables (P, T, p, etc.), we have 
five new variables associated with the convection problem: Fad, 
Fony, V, Vi, and V. However, we also have five equations: (8.1.1), 
(8.1.3), (8.1.6), (8.1.7), and (8.1.12). Thus, all the variables can 
be computed (although it takes a bit of math). The mixing length 
lm is a free parameter, and must be chosen reasonably to match 
observations as best as possible. 


Solving the Convection Equations 


The five convection equations are 
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To solve these equations, we first consolidate the constants and 
location condition variables into two dimensionless quantities. 
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and 
W = Vriad — Vad (8.1.14) 


All the quantities in these variables are known from the stellar 
structure, and can be taken as constants. 


We can then proceed to eliminate the variables. If we substitute 
the expression for velocity (8.1.6) in (8.1.12), then 
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Next, we can use (8.1.1) and (8.1.3) to eliminate Fraa 
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and use this with (8.1.7) to eliminate Funy 
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Finally, we can combine (8.1.15) and (8.1.17) by noting that if we 
re-write (8.1.15) as 
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then the equation is actually a quadratic in (V — Vive 2 with the 
(positive) solution 
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if we assign 
€=U*+V—-Vaa 
then (8.1.18) is simply 
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Now if we substitute this into (8.1.17), then 
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This is straightforward cubic equation in €, that should have only 
one real root. Since U and Vaq are both known quantities the 
solution of this equation gives V, which in turn, can be used to 
calculation Faq and Fony. 


When is Convection Important? 


The quantity U is a measure of how efficient convection is at 
transporting energy. The convective cells have two characteris- 
tic timescales: one giving the length of time required to move the 
distance lm and deliver its energy, and the other describing how 
long it takes for the cell to lose its heat via radiation. The former 
timescale can be estimated from the free-fall equation 
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The latter can be derived using our estimates of the total heat 
content (8.1.8) and rate of heat loss (8.1.9) of each blob. 
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Thus, to order unity, U gives the efficiency of energy transport. 
If U —+ 0, then the timescale for heat transport is much less 
than that for heat loss, and V should equal Vaq. This can be 




















seen via (8.1.13) and (8.1.15), since for small U, Vi © Vaa, and 
thus V ~% Vaa. This occurs in dense regions (i.e., in convective 
cores). In the outer envelopes of stars, however, the convective 
blobs can radiate a significant amount of their heat while they are 
in transit. In this case, U —> œ, so through (8.1.15), V — Vi, 
and therefore through (8.1.17), V — Vraa. The full mixing length 
theory is usually only needed in the outer parts of cool stars, where 
superadiabatic conditions can exist i.e., where both radiative and 
convective flux transportation are important. Most of the time, 
one or the other dominates, so you can use Vraa OF Vaa, Whichever 
is appropriate. 





The results you get for convective energy transport depend on the 
mixing length. In the literature, mixing length is not quoted; in- 
stead the value usually given is 
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Traditionally, œ ~ 1; currently, the “best fit” seems to be a ~ 1.7. 
(Note that this violates one of the tenants of mixing length theory.) 


Convection theory is only an approximation. Convection similar 
to that occuring in stars, cannot be studied in the laboratory. If 
we define the ratio of radiative and thermal diffusion as 
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and let v represent the viscosity of the gas, then the Rayleigh 
number 
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is a dimensionless quantity that represents how well convection 
is proceeding. For laboratory fluids, Ra ~ 10t!, but in stars, 
Ra ~ 107!. Thus, our experiments are 10 orders of magnitude 
away from reality. (In fact, the closest driven convection on earth 
is convection heating in the air above deserts.) 


Semi-Convection 


A phenomenon that takes place in the cores of massive stars is semi- 
convection. Consider a star with a convective core within which 
nuclear reactions are proceeding, and a radiative envelope. In the 
convective core, Vraq > Vaa; outside the core, Vraa < Vaa, and, 
at the border between the two regions, Vraa = Vaa. AS always, 
Vraa is given by 
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Now assume that the major source of opacity is electron scattering, 
and that the matter is fully ionized. (Since we are talking about 
the core of a massive star, these are excellent assumptions.) As 
nuclear reactions proceed in the core, the hydrogen content of the 
region will decrease. Since the electron scattering opacity for fully 
ionized matter is given by 


Ke = 0.2(1+ X) (6.1.7) 


this results in a decrease in the stellar opacity. 
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Now let’s examine the behavior of Vaq at a distance € from the 
convective core. By definition, for convection to occur 





Vradi > Waa inside the core 


and 
Vradjo < Vaq Outside the core 


Therefore, since Vaq shouldn’t change much over a distance e€, 
Vrad,i > Vrado by definition of convection 


However, Vraa is directly proportional to the hydrogen content, 
and X; < Xo. Thus, 


Vradi < Vrad,o from the opacity 


Thus there is a fundamental contradiction. 





Another way to look at the problem is to consider the temperature 
and density difference between a convective blob and its surround- 
ings. Recall that for a region to be dynamically stable, 
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where A represents the difference between the region internal to 
the blob and the blob’s surroundings. For positive displacements 
of Ar, the blob will be denser than its surroundings, and will sink 


back to its original position. 





or 











Now consider the temperature difference between the blob and the 


surrounding. 
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If the blob is in pressure equilibrium (given) and the stellar region 
is chemically homogeneous, then from (3.2.2) 
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The negative sign indicates that the inside of the blob will be cooler 
than its surroundings. It will therefore gain some energy by radi- 
ation, which will result in the blob having a slightly lower density, 
and a slightly weaker restoring force due to buoyancy. Any motion 
of the blob will eventually be damped out. 


Next, consider the situation where the surrounding medium is 
chemically inhomogeneous. We again assume that the region is 
dynamically stable, with Ap/Ar > 0. However, in this case 


=E a 


= — ell 
T dp ôu eat) 


Consequently, if 
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then blobs that have a positive Ar displacement will have AT > 
0. These blobs will be hotter than their surroundings and will 
therefore lose heat as they move. When this happens, their density 
will increase, the restoring force will increase, and their oscillation 
will slowly increase in amplitude. The blob will be vibrationally 
unstable; the result will be a slow mixing, which will eventually 
destroy the chemical gradient. This semi-convection will occur 
when 
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Note that from (8.2.1) and (8.2.2), it will only occur when the mean 
molecular weight is decreasing rapidly with radius, or equivalently, 
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This occurs in high mass stars with convective cores. In these ob- 


jects, the convective region gets smaller with time, leaving behind 
a chemical gradient. 
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